The fate of non-trivial entanglement under gravitational collapse 



Eduardo Martin-Martinez 

Institute for Quantum Computing, Department of Physics and Astronomy, and Department of Applied Mathematics, 
University of Waterloo, 200 University Avenue W, Waterloo, Ontario, N2L 3G1, Canada 



o 

^— > 

O 

O 



Oh 

> 

(N 
> 
CO 

\o 

in 
o 

CN 



Luis J. Garay 

Dcpartamento de Ffsica Teorica II, Universidad Complutense de Madrid, 28040 Madrid, Spain 
Institute de Estructura de la Materia, CSIC, Serrano 121, 28006 Madrid, Spain 

Juan Leon 

Institute de Ffsica Fundamental, CSIC, Serrano 113-B, 28006 Madrid, Spain 

Abstract. We analyse the evolution of the entanglement of a non-trivial initial quantum field state (which, for simplicity, 
has been taken to be a bipartite state made out of vacuum and the first excited state) when it undergoes a gravitational 
collapse. We carry out this analysis by generalising the tools developed to study entanglement behaviour in stationary 
scenarios and making them suitable to deal with dynamical spacetimes. We also discuss what kind of problems can be 
tackled using the formalism spelled out here as well as single out future avenues of research. 
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The question of how entanglement behaves in non-incrtial frames and in curved spacetimes has been around for already 
a fairly long time. There are many works that centre in the study of uniformly accelerated observers (among many 
others [TJ [21 El HI El El El [HI [HI HO]), or in the background of a stationary eternal black hole .11]. There are also some 
studies involving entanglement dynamics in expanding universe scenarios which have shown that the interaction with 
the gravitational field can produce entanglement between quantum field modes [HI [T3] . 

Focusing on the problem of gravitational collapse, previous works in the literature analysed the correlations between 
the outgoing and infalling modes in a gravitational collapse when the initial state is the vacuum (see for example 
[HI [TBI [TrJl ITT1 IT3"] again among many others) . 

In this work we consider the following more involved but central issue as far as the behaviour of entanglement in 
a dynamical spacetime is concerned. In the asymptotic past the field lives in a flat spacetime and its state has some 
degree of quantum entanglement between two of its modes. Then at some point, gravitational collapse occurs. The 
collapse makes the observers of the field unable to access the full state due to the formation of an event horizon. This 
has an impact on the entanglement that any observer of the field state can acknowledge. 

Studying this sort of problems is interesting from many perspectives apart from understanding how quantum 
correlations behave in dynamical curved spacetimes. Quantum entanglement plays a key role in black hole 
thermodynamics and the fate of information in the presence of horizons. Also, the study of the behaviour of non- 
trivial quantum entanglement in gravitational collapse may be useful for analog gravity proposals that aim at making 
use of this entanglement as a resource to check genuinely quantum effects derived form the formation of a horizon [17j . 
In general, this will arguably constitute a rather difficult exercise. However, inspired by tools developed to study the 
effect of accelerations on quantum entanglement, it may be possible to shed some light on this problem. 

In the study of quantum entanglement from non-inertial perspectives, i.e. in the context of relativistic quantum 
information, it was not until relatively recently that the physical meaning of the so-called 'single mode approximation' 
was analysed in detail. This approximation was introduced in [TJ [18] . It consisted in assuming that the Bogoliubov 
transformations between Minkowski and Rindler modes did not mix frequencies. In 2010 appropriate procedures to 
construct inertial modes which transform to monochromatic Rindler modes were introduced in the context of Relativistic 
Quantum Information for the accelerated scenario [TJ as well as in the stationary Schwarzschild scenario [TT]. Nowadays 
we are taking steps towards the analysis of localised states |19j and entanglement behaviour in modes contained in 
cavities [9l [10] . Although there are still a number of open questions about the analysis of localised field states and their 
possible experimental implementability, the two milestones [Tlll8| and [7] have enabled us to understand better the way 
in which entanglement behaves from non-inertial perspectives. And so, we analyse here the fundamental and qualitative 
effect of a dynamical gravitational collapse on bipartite entanglement contained in non-trivial quantum field states (that 
involve vacuum and excited states) prior to the collapse. 

In Sec. 2 we introduce the basic spacetime and quantum-field ingredients and tools to analyse the fate of 
entanglement in a gravitational collapse scenario. Section 3 is devoted to the study of the evolution of the entanglement 
of a specific bipartite quantum field state made out of vacuum and an excited state (to our knowledge, this is the first 
time that this kind of non-trivial entanglement in a dynamical spacetime is analysed) . Section 4 contains the conclusion 
and some lines of future research in this context. 

2 Gravitational Collapse 

We will consider a certain maximally entangled state of two modes of the field. This entangled state lives in a spacetime 
that is originally flat. At some point, a perturbation is produced causing the spacetime to undergo a process of 
gravitational collapse. This scenario would very well describe the process of an astrophysical stellar collapse: Prior to 
the collapse the density of a star is small enough to consider that the spacetime is approximately flat. At some point, 
the internal forces that kept the star from collapsing fail to counter the gravitational interaction and the star collapses. 
If nothing stops the collapse, it will reach a point in which an event horizon is formed. 

Let us consider the following metric written in terms of ingoing Eddington-Finkelstein coordinates as 



where r is the radial coordinate, v is the ingoing null coordinate, and M(v) = m6(v — vq). For vq < v this is nothing but 
the ingoing Eddington-Finkelstein representation for the Schwarzschild metric whereas for v < vq it is just Minkowski 
spacetime. This metric represents a radial ingoing collapsing Shockwave of radiation and it is called Vaidya metric 
(described schematically in Fig. [I]). This metric is a solution to the Einstein equations (see for instance. Ref. [15]) that, 
in spite of its simplicity, describes very well the gravitational collapse scenario that we want to analyse. Refinements 
of the model to make it more realistic only introduce subleading corrections. In particular this model captures, up to 
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Figure 1. Carter-Penrose diagrams for gravitational collapse: Stellar collapse (solid line) and ingoing radiation shock- 
wave (dashed line) giving rise to Vaidya spacetime. 



subleading corrections, the more realistic collapse of a matter cloud. Let v n = vq — 4m be the coordinate of the last null 
ray that escapes to the future null infinity y + and hence that will eventually form the event horizon (see Fig. §. 

Consider now a state of a scalar quantum field. We need to introduce convenient bases of solutions to the Klein- 
Gordon equation determined by their behaviour in the different relevant regions of this collapsing spacetime. For this, 
we will follow a standard procedure (see e.g. [T5]V 

We first define the 'in' basis of ingoing positive frequency modes, associated with the time parameter v at the null 
past infinity J^ - : 

u™ ~ (4 7 rr v ^)- 1 e- i ^, at J'. (2) 

Second, we define another basis in a Cauchy surface in the future. In this case, the asymptotic future J^ 4 " is not a 
Cauchy surface in itself, so we need to consider also the future event horizon Let us begin with the 'out' modes 

defined as being outgoing positive- frequency in terms of the natural time parameter ?7 out at which are 

uZ ut ~ (iTiry^)-^-' 1 ^^, at S+, (3) 

where ?7 ou t = v — 2r* ut and r* ut is the radial tortoise coordinate in the Schwarzschild region. At early times, these modes 
it° ut concentrate near v H at and behave in the following way: 

uZ ut ps (4^r v ^r 1 e- iu ( , '»- 4 ' nln ^)%H - «), at J~ , (4) 

having support only in the region v < u B , since only the rays of light that depart from v < v H will reach the asymptotic 
region .y + . The rest will fall down into the horizon. 

Finally, we use an analytical continuation argument to define the 'hor' modes at J / t? + : These will be modes that 
behave as ii° ut in the asymptotic past J?~ , but for v > v H . In other words, we define them as modes that leave the 
asymptotic past to fall into the horizon, never reaching the asymptotic future. Near the Cauchy surface these 
modes behave as 

ul m - (4^r v / ^)" 1 e +i "( , '»~ 4mln ^)c?(i> - v H ), at J~ . (5) 

By expanding the field in terms of the two sets of modes ('in' on the one hand and 'hor-out' on the other) we can 
relate the two sets of solutions via the corresponding Bogoliubov coefficients. We will provide more details below, but 
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for now let us refer to the extensive literature on this topic (see for instance [TS] ) and directly give the expression of the 
bosonic annihilation field operators in the asymptotic past in terms of the creation and annihilation operators of 'out' 
and 'hor' modes which, in the notation of [13], read: 
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where tanh r u 



The values of <p and a uu » will be given below. The 'in' vacuum, given by a™ |0) in = for all 
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positive frequencies u>, can be readily rewritten in the 'out-hor' basis as 
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where \nj) denotes a mode with occupation number n and frequency uj. 

The next task is to write the one-particle state in the past |lw); n as a linear combination of the 'out-hor' basis 
modes. If we have a monochromatic excitation in the asymptotic past, equation ([6]) tells us that it will become a highly 
non-monochromatic linear combination of 'hor' and 'out' modes. The standard well-known procedure (which we briefly 
summarise in what follows, see e.g. [T51 IT]) is to construct another basis of 'in' modes with positive-frequency in the 
past such that its Bogoliubov transformation into 'hor' and 'out' modes is diagonal in frequencies. Let us call those 
modes Uq and Uq in order to keep the notation of Ref. [7J. Note that, for all the reasons discussed above, these modes 
are intrinsically non-monochromatic in the asymptotic past, i.e. in the basis of modes it™, and that f2 is labelling the 
frequency of such modes Mq and Uq with respect to the time in the asymptotic future region 'out'. 

As suggested by the fact that |0) in is a two-mode squeezed vacuum of 'out' and 'hor' modes, let us define new 
positive-norm 'R-L' modes by the following diagonal Bogoliubov transformation from 'out-hor' basis with the form of a 
two-mode squeezing operation: 



coshMu^f 
coshM^ 01 ' 



sinhM^* 
sinh(r n K ut * 



(8) 



Taking Klein-Gordon inner products with the 'in' modes we find the following form of Uq and Uq in terms of the u" 
modes: 



[cosh(r n )a un + smh(rn)£ w n] VJdu;, 
[cosh(m)7u,n + sinh(r f2 )/3 tlJ n]*'u™dw, 
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Pun = - tanh? 



<Ln = - tanh 



(9) 
(10) 

) are the following Bogoliubov 

(11) 

r u a* n . (12) 



Therefore we see that the 'R-L' modes are purely positive frequency linear combinations of the 'in' modes and that they 
also form a complete set of solutions of the field equation in the asymptotic past. 

This relationship between the modes directly translates into a relation between the particle operators associated 
with them 



1q = cosh ro a™* - sinh r fi 



hort 



cosh rn a 



hor 



sinh rn a! 



outt 



(13) 
(14) 



Obviously, these operators annihilate the 'in' vacuum: Oq |0) in = |0) iE 



0. 



To summarise, these modes have the following features: 

(i) They share the same vacuum |0) in as the monochromatic modes a™. 

(ii) They form a complete basis of solutions to the field equation which are positive-frequency in the asymptotic past. 

(iii) They translate into a single frequency mode when expressed in the future basis. 
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Properties (i) and (ii) allow us to decompose any physical state as some combination of these modes, which makes them 
worth studying as an intermediate stage of more general cases. The third feature (diagonal Bogoliubov transformations) 
greatly simplifies the formalism, enabling us to use all the artillery already deployed in other simpler scenarios also in 
the case of stellar collapse, providing us with a nice and clear interpretation to the analysis of the entanglement in the 
asymptotic future. 

Repeating an analogous reasoning as in [TJ, we can still introduce a more general annihilation operator with this 
properties, which will be a linear combination of the two annihilation operators Oq and o,q defined above: 

Ca 9R = fe^n + ?LGtfi, (15) 

where qi, and (?r are real parameters satisfying = yJl — qk and < q^_ < 1. 

3 Entanglement behaviour 

Let us consider the following maximally entangled bipartite state in the asymptotic past, "prepared" long before collapse 
starts 

l*> = ^(|0)^|0)f n + |l)^|ln)f n ), (16) 



where the excited modes for Bob are chosen to be those generated by (15), namely, 



whereas Alice's mode can be chosen arbitrarily (that is why it is not labelled with fl). This initial state (16) will be 
observed by two observers Alice and Bob. While we will consider that Alice has unrestricted access to her partial state, 
we will assume that Bob does not because at some point the process of gravitational collapse will generate an event 
horizon, preventing him from accessing the full state. 

Although perhaps the most natural scenario would be that in which both subsystems are in the proximities of a 
stellar collapse and, hence, both observers would undergo similar processes, let us consider for simplicity that only one 
of the subsystems is going to be affected by the stellar collapse. One can think of Alice's state prepared such that it is 
a localized state living far away from the collapsing star. Alternatively we could consider that Alice can measure her 
subsystem prior to the formation of the horizon so that it cannot hinder her ability to obtain information about her 
partial state. On the other hand, Bob's knowledge about his subsystem is going to be limited because between the time 
when the state was created and the time in which he will be able to measure it, an event horizon appears, preventing 
him from accessing the full state in the future. This scenario is not devoid of physical interest. As we will show later 
on, it will allow us to focus on questions regarding quantum correlations between modes in the past and modes falling 
into the horizon. 

In these circumstances, Alice measures in the 'in' basis whereas Bob measures in the 'out' basis, having lost all 
the information contained in the modes 'hor' that are bound to fall into the forming black hole. Now, to describe the 
effective state to which Alice and Bob can access requires that we trace out the modes that become causally disconnected 
from Bob due to the formation of the horizon, i.e. the 'hor' modes: 

PA-out = tr hor (|*)(*|). (18) 

This state is non-separable and we can compute its negativity |20j as a convenient quantifier of quantum entanglement. In 
simple words, when we compute the negativity of the density matrix pA-out we obtain a measurement of the correlations 
between Alice's fully accessible state and the modes that reach the asymptotic future escaping the collapse. 

It is legit then to ask about what would be the quantum entanglement between Alice's state and those modes that 
will not make it to the asymptotic future because they will fall into the incipient horizon, becoming trapped into the 
black hole. The state whose separability we would have to analyse would then be 

pA-hor = tro» t (|*>(*|). (19) 

The procedure to compute the negativity is the following: first we express |^) in the basis of 'hor-out' modes for 
Bob, which yields 



w =7i 



cosh rn 



j^- tanh " r ^ \°)t Mout l«n)iL 

' (20) 



-— ^— ]T v / n+Ttanh" r n |1>£ ( te |(n + l)n) B ut K)L + 9l |nn>* t |(n + l) n )£j 
cosh r n 
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Then the density matrices pA-har and pA-out are obtained after a simple but lengthy algebra exercise by tracing 
out the 'out' and 'hor' modes respectively from |$), 



PA-out 



£(" + (<& |1>£ l(» + l)n)L <(" + ^ll + it |1>£ Ml <1|£ (nnll) 
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«R |0)r„ Ml (l|r„ <(" + l)nll + IC l(" + l)n)l Ml 



B 

out 



|l)f n Ml (0|£<(n + l) n |* 



(21) 



where T r = tanhr^ and C r = coshrji. The expression for pA— hoi is obtained by exchanging and qi, and 'out' by 
'hor' in the equation above. 

The entanglement monotone that we will compute, the negativity |20j . is the sum of the negative eigenvalues 
of the partial transposed density matrix of the quantum state for which we want to evaluate its degree of distillable 
entanglement. To compute it we first take partial transposes in pA-out and pa-Uot (which is the transpose only with 
respect to Alice's indices). This yields 

rp2n 

— — In 

2a 2 



PA-out — 



+ Z)( n + ^ l^ta K n + (lite <(« + l)nl^t + ^ Il>ta l"n>^ t <1|^ <nn|^ t ) 



n=0 



rp2n+l 

£ ^(n+l)(n + 2)-^- to g L (|1>£ |n n )l <1|£ ((n + 2)„|^ t + |1>£ |(n + 2) n >l <1|£ MlJ 
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n=0 
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A Ml (0|f n ((n + l)o|l + |0>£ |(n + l)n)l <1|£ Ml 



- 7l ( |l)f n |(n + l)n>l (Q|£ (noil + I0)f n Ml ((n + l) n |* 



(22) 



where, as above, the expression for Pa_^ ot is obtained by exchanging q^ and qi, and 'out' by 'hor'. 

As it is the case of the accelerated-observer scenario, the diagonalisation of the infinite-dimensional partial 
transposed density matrices p^_ hor and p^_ hor can be carried out only numerically since, with the exception of the case 
<7r = 1, no block-diagonalisation can be performed. 

Figure [2] shows the result of the calculations. The negativity 'A-out' as a function of the mass of the forming black 
hole and the frequency of the probed 'out' mode is shown as solid blue lines. The negativity 'A-hor' is plotted in the 
same figure as dashed red lines. 

We see that varying the parameter q^ we are basically controlling whether Alice's state will have more quantum 
correlations with the modes that will reach the asymptotic future or the modes that will fall into the horizon. This is 
best seen in the large black hole mass limit, where there is an exact trade-off between the correlations that Alice's mode 
has with infalling and outgoing modes. 

We can also see that for very small black holes, all entanglement is completely degraded: when m no 
entanglement survives either with the modes that fall into the horizon or the modes that will reach the future. As the 
mass of the black hole (or the frequency of the probed mode) is increased, the correlations quickly become insensitive to 
the presence of the horizon, as one can expect taking a look at the quantum effects induced by gravity in the presence 
of a black hole: They become stronger as the mass of the black hole is closer to zero. This suggests that this can be 
understood in a pictorial way as a limit in which the Hawking-like radiation spoils all correlations contained in the state. 

As it is well known [21 , if Bob can only measure modes in the asymptotic future, he will see the vacuum state |0) in 
as a thermal state. Indeed, if we compute how the 'in' vacuum is seen by observers in the asymptotic future we obtain 
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Figure 2. Negativity of PA— out (solid blue line) and P/i_hor ( re( i dashed line) as a function of the product of mass of 
the black hole and probed frequency for various choices of the past modes | |15| l . Note that no entanglement survives the 
singular black hole limit (m — > 0) where the Hawking temperature is divergent. 
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Pout.i 



(coshr^) 2 f-^> 

v 7 n— 



(taring) 2n K) Jn u \ 



(23) 



This is a thermal radiation state whose temperature is Th = (Sirm) 1 . So if that happens with the vacuum state it would 
be a reasonable hand-waving argument that, if instead of the vacuum, we consider a pre-existing non-trivial entangled 



state such as that of equation (16), the thermal-like noise could impair the ability of the observers of acknowledging 



quantum correlations in the system. 

However one has to be very careful when thinking to what extent this behaviour can be naively associated to the 
Hawking thermal noise. To begin with, we are not considering the vacuum state, but rather an entangled state of field 
excitations. The process of change of basis and tracing out of the modes that fall into the event horizon is not as trivial 
as for the vacuum case. In fact, it has been shown that in the Rindler scenario |22j and beyond the so-called single mode 
approximation with some choices of the state, the accessible entanglement for an accelerated observer may behave in a 
non-monotonic way, as opposed to the first results reported in [51 [3J [7] . This is due to inaccessible correlations in the 
initial states becoming accessible to the accelerated observer when his proper Fock basis changes as acceleration varies. 



While this phenomenon was highlighted in [22] for the Rindler case, for an analogous choice of the modes (15), a similar 
behaviour would be expected in the dynamic scenario analysed here. 

Let us conclude with a note of warning: The modes analysed here as tools to study entanglement behaviour in 
gravitational collapse have very nice properties but due to their highly non-monochromatic nature and non-localisation 
they are modes that can arguably be difficult to prepare and measure in an hypothetical experiment. This said, this 
tool will allow us to simplify the calculations so that we can extract fundamental results in settings in which other 
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techniques have proven not operational, much in the same fashion as the introduction of Unruh modes [7J has allowed 
progress in our understanding of quantum correlations from non-inertial perspectives. One has to keep in mind that the 
modes used here share some fundamental properties with the standard monochromatic 'in' modes, and that they form 
a complete basis of solutions to the field equations. This means that any physically conceivable state can be expressed 
as a superposition of modes as the ones studied here. 

4 Conclusion and Future research 

We have analysed the behaviour of quantum entanglement present in some initial field state when a gravitational collapse 
occurs in the background. The quantum correlations in that state are perturbed by the formation of the event horizon, 
mixing the quantum state originally prepared, and therefore, degrading the original correlations. We have done so by 
adapting the tools developed in the analysis of entanglement in the context of the Unruh-Hawking effect [3J [IT] to go 
beyond what was known as 'single mode approximation' [7J. 

We have shown, that similarly to the infinite acceleration limit in accelerated scenarios [2J O [7] and similar to the 
stationary eternal blackhole scenario entanglement is completely degraded when we consider singular black holes 
(to — > 0) for which the Hawking temperature diverges. 

A trivial extension of the results obtained here is considering that the appearance of an event horizon affects the 
ability to access the full state for both Alice and Bob. In these cases, and for maximally entangled states of a scalar 
field, entanglement will be arguably degraded more quickly than in the case where only Bob is affected by the collapse, 
much in a similar way as it happens in the acceleration scenario |23j . Extending this result to the case where both 
observers measure after the horizon is created is somewhat straightforward with the tools developed here, being mainly 
a matter of a more complicated calculation, whereas the results are arguably going to be qualitatively the same. 

The next natural step is to introduce localised measurements that will endow the entanglement degradation 
phenomena reported here with operational meaning. Using for example, localised projective measurements as in |19j 
more physical scenarios can be analysed. 
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